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Preface 


This  study  of  the  atomic  motion  in  a  standing  wave  was 
motivated  by  the  challenge  presented  in  understanding  the 
basic  physical  processes  involved.  The  study  of  resonance 
radiation  pressure#  while  dating  back  to  the  origins  of 
quantum  mechanics#  is  still  a  new  and  unexplored  field.  I 
am  deeply  indebted  to  my  advisor#  Major  Richard  J.  Cook#  f 
sparking  my  interest  in  this  field#  and  for  helping  me  to 
understand  the  nature  of  the  investigation.  I  would  also 
like  to  express  my  gratitude  to  all  my  friends  in  Dayton# 
Ohio,  who  helped  me  through  the  difficult  times  with  their 
support  and  prayers. 
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Abstract 


The  motion  of  an  atom  in  an  intense  resonant  standing 
laser  wave  is  studied  using  the  Ehrenfest-Bloch  equations. 

Two  approaches  are  taken  to  determine  the  radiation  or  light 
pressure  force  on  the  atom.  The  first/  an  attempt  to  find  a 
Fourier  series  solution  to  the  optical  Bloch  equations 
failed,  but  Fourier  coefficient  recursion  formulas  were 
found.  The  second,  a  numerical  calculation  of  the  Bloch 
variables  and  the  radiation  force  was  successful  in 
obtaining  detailed  plots  of  the  instantaneous  values.  The 
effect  of  assuming  constant  atomic  velocity  is  studied,  and 
Fast  Fourier  Transform  descrete  harmonic  frequency  spectra 
of  one  of  the  Bloch  variables  and  the  light  pressure  force 
are  obtained.  Finally,  the  Ehrenfest-Bloch  equations  are 
compared  with  other  theories  of  resonant  radiation  pressure. 


ATOMIC  MOTION  IN  A  STANDING  WAVE 


I .  I n  troduc  tion 

With  the  advent  of  hi<jh  power  tunable  laser  systems, 
the  study  of  the  effect  of  radiation  pressure  on  atoms  r.ear 
atomic  resonant  frequencies  has  developed  into  a  fertile 
field  of  investigation.  of  particular  interest  in  recent 
years  is  the  study  of  the  motion  of  an  atom  in  an  intense 
standing  laser  wave.  The  light  pressure  of  such  standing 
waves  may  provide  a  means  to  damp  out  unwanted  thermal 
motion  in  a  beam  of  particles,  or  to  cool  and  trap  single 
neutral  atoms  within  traps  of  laser  light.  With  such 
setups,  doppler-free  spectroscopy  could  be  performed  on 
single  trapped  atoms.  Single  atom  experimentation  may 
become  the  norm  of  modern  experimental  physics.  A  trapped 
neutral  atom  could  also  be  used  as  an  oscillator  to  create  a 
highly  accurate  atomic  time  base  for  physical  measurements. 
The  potential  uses  of  light  pressure  to  change  the  motion  of 
an  atom  is  a  long  and  very  imaginative  list. 

But,  to  be  able  to  use  the  effects  of  radiation  force 
to  perform  these  investigations,  it  is  necessary  to 
characterize  and  understand  what  happens  between  the  light 
wave  and  the  atom.  Several  theorems  of  light  pressure  have 
been  developed  which  describe  interaction  between  the  light 
wave  and  the  atom.  This  thesis  describes  the  use  of  the 
Ehrenf est-Bloch  equations  to  characterize  the  atomic  motion 


in  an  intense  standing  laser  wave. 

In  this  work/  an  attempt  was  made  to  analytically  solve 
the  optical  Bloch  equations  using  a  Fourier  series  solution. 
Although  it  was  unsuccessful  in  obtaining  a  complete 
solution  to  the  Bloch  equations,  a  set  of  Fourier 
coefficient  recursion  formulas  were  found  that  can  be  used 
to  study  the  coupling  between  terms  in  the  series. 

In  addition  to  the  analytical  efforts,  the  Ehrenfest- 
Bloch  equations  were  solved  numerically  to  get  a  picture  of 
the  atomic  motion  as  a  function  of  time  and  field  intensity. 
Detailed  plots  were  created  that  describe  the  internal 
motion  of  the  atom  with  the  Bloch  vector,  and  the  external 
motion  under  the  influence  of  the  light  pressure  force.  As 
a  part  of  the  study,  the  assumption  of  constant  atomic 
velocity  and  its  effect  on  the  results  of  the  numerical 
calculations  was  investiga ted .  Also,  a  Fast  Fourier 
Transform  was  performed  on  the  Bloch  variable  ZL  ,  and  the 
light  pressure  force  to  determine  their  dependence  on  higher 
harmonic  frequencies  as  a  function  of  field  intensity. 

The  results  obtained  from  using  the  Ehrenfest-Bloch 
equations  are  comparable  to  the  published  results  obtained 
through  other  theories  of  the  light  pressure  forces.  But, 
by  using  the  Ehrenfest-Bloch  equations,  one  obtains  a 
detailed  instantaneous  view  of  the  atom  as  it  moves  through 
the  standing  wave.  This  view  is  a  significant  step  toward 
the  complete  understanding  of  the  interaction  between  light 
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Radiation  Force 

Laser  light  tuned  to  a  frequency  near  the  resonant 
transition  frequency  of  an  atom  is  strongly  scattered  by  the 
atom.  In  the  scattering  process,  momentum  is  transfered 
from  the  field  of  the  light  to  the  atom.  The  rate  of  this 
momentum  transfer  is  called  radiation  force  or  light 
pressure  on  the  atom.  Light  pressure  is  a  composition  of 
two  forces  on  the  atom:  the  spontaneous  or  longitudinal 
force,  and  the  dipole  or  gradient  force. 

The  spontaneous  force  results  from  the  absorption  of 
the  light  wave  which  excites  the  atom,  and  the  subsequent 
decay  of  the  atom  to  the  ground  state.  Because  the 
spontaneous  emission  occurs  in  a  completely  random 
direction,  the  atom  feels  no  net  change  in  momentum  due  to 
the  emission.  The  momentum  change  is  due  to  the  absorption 
of  the  light  wave  and  thus  occurs  in  the  direction  of  light 
propagation,  hence  the  name,  longitudinal  force. 

The  dipole  force,  on  the  other  hand,  is  a  result  of  the 
interaction  between  the  induced  atomic  dipole  moment  and  the 
amplitude  gradient  of  the  electric  field  of  the  light  wave. 
This  force  is  characterized  by  its  direction  toward  or  awa 
from  increasing  field  strength  depending  on  whether  the 
laser  is  tuned  to  a  frequency  below  or  above  the  resonant 
transition  frequency  of  the  atom.  When  the  laser  is  tuned 
exactly  on  resonance,  the  atom  experiences  no  dipole  force. 
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Because  the  of  the  profile  of  a  laser  beam,  this  force  is 
usually  directed  perpendicular  to  the  direction  of  light 


propagation.  Thus,  the  dipole  force  is  often  called  the 
transverse  force. 

Atomic  Cooling  in  Resonant  Radiation 

The  study  of  atomic  motion  in  resonant  or  near-resonant 
radiation  has  received  considerable  interest  in  recent  years 
since  the  proposal  by  T.W.  Hansch  and  A.L.  Schawlow  to  use 
light  pressure  force  to  cool  atoms  by  placing  them  within  a 
standing  light  wave  with  a  frequency  tuned  slightly  below 
the  resonant  transition  frequency  of  an  atom  [6].  In  their 
proposal,  the  thermal  motion  of  the  atoms  will  cause  the 
frequency  of  the  light  to  be  doppler  shifted  such  that  one 
of  the  two  coun ter- pcopaga ting  waves  of  the  standing  wave 
closer  to  resonance,  while  the  other  would  be  moved  farther 
from  resonance.  This  would  cause  the  atoms  to  absorb  more 
momentum  from  the  wave  closer  to  resonance,  and  thus 
experience  the  spontaneous  force  opposing  their  motion.  In 
this  way,  the  motion  of  the  atoms  in  the  standing  wave  would 
be  damped  out  and  the  atoms  cooled. 

This  relatively  simple  idea  of  applying  light  pressure 
to  cool  atoms  led  to  the  development  of  detailed  theories  of 
the  light  pressure  force  [3],  [5],  [7],  [10],  [13].  These 
theories  characterized  the  nature  of  the  light  pressure 
forces  and  found  them  to  be  more  complex  than  first 
envisioned  by  Hansch  and  Schawlow. 


Neutral  Atom  Trapping 

As  the  theories  of  light  pressure  were  being  developed/ 
it  became  apparent  that  it  was  theoretically  possible  to  use 
radiation  pressure  to  create  a  potential  well  into  which  an 
atom  could  be  placed.  The  idea  of  a  "light  trap"  to  hold  an 
extremely  cold  atom  within  a  beam  of  light  spurred  interest 
in  the  study  of  the  motion  of  an  atom  in  an  intense  standing 
laser  wave  [2],  [9].  As  discussed  in  the  introduction,  the 
successful  cooling  and  trapping  of  a  neutral  atom  would 
revolutionize  experimental  physics,  and  indeed,  is  one  of 
the  most  ambitious  goals  of  modern  physics  today. 


Doppleron  Resonances 

In  the  development  one  theory  of  the  light  pressure  in 
a  standing  wave,  E.  Kyrola  and  S.  Stenholm  predicted  the 
existence  of  multiphoton  or  Doppleron  resonances  to  occur  at 
when  the  atomic  velocity  was 


=.  ± 


k(2n-H) 


This  was  described  as  the  result  of  absorbing  n  +• I  photons 


from  one  of  the  two  counter-propagating  light  waves  that 
make  up  the  standing  wave,  and  emitting  n  photons  into  the 
other.  The  name,  Doppleron,  was  applied  because  the  process 


is  caused  by  the  Doppler  shift  of  the  field  for  the  moving 
atom.  These  resonances  were  found  to  become  significant  for 


large  field  strengths. 


Calculations  from  a  theory  of  resonance  light  pressure 
by  V.G.  Minogin  and  O.T.  Serimaa  showed  that  the  Doppleron 
resonances  predicted  by  Kyrola  and  Stenholm  appeared  in  the 
form  of  variation  in  the  light  pressure  force  [11]. 

They  presented  the  light  pressure  force  as  the  sum  of 
two  orthogonal  components:  the  longitudinal  force  F„  ,  and 
the  transverse  or  gradient  force  FA  .  These  components 
were  expressed  in  the  form  of  Fourier  expansions: 

F„  =  F,*  +  £  ( F„*n  c°5  2nkx  +  F*nsin2nkx)  (2a) 

Fx  =  F*  +  £  (  FaCb cos  2nkx  +  Fx‘nsinZnkx)  (2b) 

where  the  coefficients  are  a  rather  complicated  set  of 
convergent  partial  fractions.  For  very  high  intensities  and 
and  low  atomic  velocities,  the  force  actually  changes  sign. 
Instead  of  being  a  damping  or  cooling  force,  the  light 
pressure  force  becomes  an  accelerating  or  heating  force. 

The  force  resonances  found  by  Minogin  and  Serimaa 
are  centered  approximately  at  the  velocities  predicted  in 
equation  (1),  shifted  toward  greater  velocities  and  power 
broadened  as  the  intensity  of  the  radiation  field  increases. 

Figure  (1)  is  a  plot  of  the  time  averaged  longitudinal 
force  calculated  from  the  expansion  of  the  Fourier 


coef f icients. 


Figure  1.  Force  as  a  function  of  atomic  velocity 
for  various  field  intensities  [11], 


Ehrenf est-B loch  Equations 

The  difficulty  of  working  with  the  statement  of  the 
light  pressure  force  developed  by  Minogin  and  Serimaa  lies 
in  the  complexity  of  the  partial  fraction  solution  for  the 
Fourier  coefficients  of  the  force. 

A  simplified  theory  of  near  resonance  light  pressure 


was  developed  by  R.J.  Cook  [3J.  This  theory,  based  on 
Ehrenfest’s  Theorem  and  the  optical  Bloch  equations  provides 


a  simple  statement  of  the  radiation  force  through  a  system 
of  four  coupled  equations  called  the  Ehrenfest-Bloch 
equations 


u  =  -  +  ait 

XT  -  -  &ru  -  4Av-  +  £1  XT 

=  —  Slv  -  A (o«r  +  \) 

and 

p  = 

Equations  (3)  are  the  optical  Bloch  equations  for  a  two- 
level  atom  in  the  Rotating  Wave  Approximation  (Appendix  A), 
and  Equation  (4)  is  the  equation  of  motion  of  the  atom 
obtained  from  Ehrenfest' s  Theorem. 

For  the  case  of  a  standing  wave,  the  electric  field  is 
a  function  of  position,  £00  =  £a  Cjo*  kx  ,  and  the  phase  is 
given  by  6  —  0  .  Therefore,  the  Rabi  frequency  becomes 
SI.'**  Si. ,  t»s  lo  where  SL»  —  ^*-£» / .  The  radiation  force 

can  then  be  written  as 

P  =  -i'fck-fU'W  Sin  kx  (5) 

If  the  velocity,  V  ,  of  the  atom  remains  constant  (a 
reasonable  approximation),  then  ><  =■  and  t^e  Ehrenfest- 
Bloch  equations  become 


(3a) 

(3b) 

(3c) 


'Ll  -  -  ikv.  +  azt 

ir  -  —  &V.  -  fAv  +  il.coskvt  yr 


(6a) 


(6b) 

IaT  —  — _n#GOS  kvt  IT  -  A  (w  +0  (6c) 

and 

F  =  -  i-Fik-0..  U  sin  kvt  (7) 

The  light  pressure  force  is  determined  by  solving  the  Bloch 
equations  and  substituting  the  value  of  U  into  equation 
(7). 

The  Weak  Field  Case 

The  Ehrenf est-Bloch  equations  can  be  solved 
analytically  for  the  case  of  the  weak  field  excitation  [11] 
In  this  case,  the  atom  will  remain  near  the  ground  state 
(W«-l)  and  the  Bloch  equations  can  be  written  as 


lA  -  -  ^ A'U  4-  A LT  (8a) 

xr  =■  -  A?X  -  jf\ir  -  i7.0cos  kvt:  (8b) 

XJ '  —  0  (8c) 

The  solution  for  this  set  of  equations  was  found  to  be 

IX  =  OX  M>s  kvt  4-  fS  sin  kvt:  ( 9a ) 

t  cos  kvt  +  8  sin  kvt  :  -  ] 


%r  - 


(  Qh  \ 


where  the  coefficients  can  be  found  by  solving 


*  _  ±M  -  A Jb 


(3  = 


kv 

-■kA«-  +  aK 


kv 


„  _  M  +  xAS 
* - k7 - 

c  -  -a°<.--z Ay  -  -Q.« 

6  kv 


( 10a ) 


(10b) 


(10c) 


(10d) 


The  force  is  given  by  substituting  the  solution  for  U.  into 
equation  (7)  to  get 


F  =  - iSl^k  (<*  cos kvt.  +  /S  sin  kvt)  sin  kvt 


(ii) 


Or,  when  the  force  is  averaged  over  time. 


F  =  '^Sl.*kf3  (12> 

where 

^  =  _ -  a  ilpAkv _ 

[  a1  -  (kv)1]1  +  ttA)1  [  (if\)h  2&  *-2(kv)1]  U3) 

This  shows  that  in  the  weak  fieLd  case,  the  direction  of  the 
force  depends  on  the  product,  AV  .  A  negative  force 
opposes  the  motion  of  the  atom  and  will  slow  it  down,  while 
a  positive  one  is  in  the  direction  for  the  atomic  motion  and 
will  speed  it  up.  This  is  in  agreement  with  the  weak  field 
results  by  Minogin  and  Serimaa. 
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The  first  part  of  the  thesis  research  was  spent 
analytically  manipulating  the  optical  Bloch  equations  to  see 
if  they  could  be  solved  with  a  Fourier  series.  This  series 
was  then  to  be  investigated  to  determine  the  dependence  of 
the  time-average  1  force  on  the  magnitude  of  the  electric 
field  of  the  standing  laser  wave. 

This  effort  was  unsuccessful  in  obtaining  a  Fourier 
series  solution  to  the  Bloch  equations.  Several  recursion 
formulas,  however,  were  found  for  the  coefficients  of  the 
Fourier  series.  This  section  describes  the  Fourier 
coefficients  and  their  use  in  approximating  the  weak  field 
solution . 


Fourier  Series  Solution 

The  weak  field  solution  was  found  by  approximating  the 
optical  Bloch  equations  with  the  assumption  that  the  atom 
remained  near  the  ground  state.  While  this  may  be  true,  no 
information  about  the  strength  of  such  a  weak  field  can  be 
determined  from  the  initial  assumptions.  Using  the  Fourier 
solution  to  solve  the  complete  Bloch  equations,  the  weak 
field  case  can  be  investigated  by  studying  this  solution 

From  equations  (9),  the  complete  solution  of  the  Bloch 
equations  could  be  assumed  to  be  a  Fourier  series  in  which 
the  higher  order  coefficients  become  negligible  in  the  weak 
field  case.  Such  a  solution  would  be 
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where  the  coefficients  would  have  to  be  determined 
the  complete  solution.  By  substituting  the  assumed 
into  equations  (6),  one  obtains 
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to  obtain 
solution 
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Using  the  trigonometric  identities, 
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equation  (20)  becomes 
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The  summation  indices  in  the  last  two  terms  in  equat 
can  be  changed  to  read 
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Equation  (22)  thus  becomes 
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And  by,  equating  coefficients, 
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Similarly/  equation  (15c)  can  be  rewritten  as 
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And  from  equating  coefficients/ 
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Equations  (17)  -  (19),  (28)  -  (32),  and  (34)  -  (38)  form  a 
coupled  system  of  recursion  formulas  for  determining  the 
Fourier  coefficients  for  the  complete  solution  to  the  Bloch 
equations.  While  they  were  not  solved  for  each  coefficient 
during  this  study,  they  are  still  useful  to  provide  some 
insight  into  the  nature  of  the  Fourier  solution. 
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The  Truncated  Fourier  Series 

For  the  weak  field  case,  the  Bloch  equations  were 
approximated  according  to  the  assumption  that  the  atom  would 
remain  in  or  near  the  ground  state.  Equations  (9)  list  the 
analytical  solutions  in  the  weak  field  case. 

If  the  recursion  formulas  obtained  for  the  complete 
Bloch  equations  depict  the  coefficients  of  the  Fourier 
series  solution,  then  they  would  have  to  be  the  same  as  (9) 
in  the  limit  of  a  weak  field.  Assuming  that  the 
coefficients  of  the  higher  harmonic  terms  are  negligible  in 


the  weak  field  case,  the  coefficients  can  be  written 


IV 


Numerical  Solution 


The  second  area  of  effort  in  this  thesis  was  to  use  the 
Ehrenf est-Bloch  equations  to  determine  the  instantaneous 
motion  of  the  atom  as  it  moves  through  the  standing  wave. 

The  computer  programs  that  were  used  were  written  in  Fortran 
77  and  can  be  found  in  Appendix  B. 

Because  of  the  way  the  Ehrenf est-Bloch  equations  were 
derived  for  the  case  of  a  strong  standing  wave,  two  methods 
were  used  to  program  the  equations.  The  first  method  kept 
the  atomic  velocity  constant  as  in  equations  (6)  and  (7). 

The  second  method  allowed  the  velocity  to  vary  as  it  would 
in  an  actual  experiment.  Both  methods  provided  useful 
information  about  the  atomic  motion  in  a  standing  wave. 

Constant  Velocity  Program 

The  Ehrenf est-Bloch  equations  were  first  programmed 
with  the  assumption  that  the  atomic  velocity  was  constant. 

The  values  used  for  the  Einstein  A-Coef f icient  and  the 
atomic  transition  frequency  were  A=  10*  s&c  1  and 

io“*  sec"'  respectively.  These  values  were  chosen 
arbitrarily  and  are  not  actual  values  for  real  atoms.  The 
detuning  was  chosen  to  place  the  laser  frequency  one 
A-coe f f ic ien t  below  the  atomic  transition  frequency. 

The  Bloch  equations  were  programmed  exactly  as  they 
appear  in  equation  (6).  The  derivatives  of  the  variables 
were  calculated  at  the  beginning  of  each  time  step,  then  the 


time  (x  IO'ssec) 

Plots  of:  (a)  field  strength, 

(b) ,  (c) ,  and  (d)  Bloch  variables, 

and  (e)  radiation  force  (dynes)  for 
constant  velocity  v  =  6*  IO*  cm/set  • 


time  was  incremented  and  the  new  values  for  the  variables 
were  found  by  a  linear  calculation  across  the  time  step. 

The  size  of  the  time  increment  that  was  used  was  )0  "  sec. 

Figures  (2)  and  (3)  are  pi  s  of  the  Bloch  variables 
and  the  radiation  force  that  were  calculated  for  .fio®  5A 
with  the  constant  ve  loci  ties  ,  v  =  2*io5  w/i«c  and  6  *  io*  on/sec 
respectively.  From  these  plots,  it  can  be  seen  that  the 
largest  contribution  to  the  light  pressure  occurs  before  and 
after  the  atom  passes  through  a  node  in  the  standing  wave. 
This  is  because  the  dipole  force  depends  on  the  gradient  of 
the  field  amplitude  which  is  the  greatest  at  the  nodes. 

Close  observation  of  Figures  (2e)  and  (3e)  show  that  the 
light  pressure  force  for  these  velocities  and  field  strength 
will  not  average  to  zero.  In  fact,  the  average  force  from 
Figure  ( 2e )  is  a  small  positive  quantity,  while  the  average 
force  from  (3e)  is  a  small  negative  quantity.’ 

It  can  also  be  noted  that  the  frequency  of  V.  and  V 
is  equal  to  the  frequency  of  the  standing  wave  as  it  is 
experienced  by  the  moving  atom,  or  the  quantity,  kv  .  The 
force  and  the  inversion  seem  to  have  frequencies  of  twice 
this  frequency  or  Zkv  .  And  finally,  Fcabi  oscillations  are 
observable  in  Figure  (3)  for  V  and  "W"  as  a  ringing  in 
the  Bloch  variables. 

Effect  of  fx  on  the  Radiation  Force 

In  equation  (7),  the  magnitude  of  the  light  pressure 


force  for  ttie  case  of  an  atom  moving  through  a  standing  wave 
depends  on  the  amplitude  of  the  standing  wave  explicitly 


time  (  x  10"*  sec.) 


Figure  4.  Plot  of  XI  for  various 
field  strengths. 


through  the  Rabi  flopping  frequency,  £1.  = ,  and 
implicitly  through  the  Bloch  variable  U  . 

Figure  (4)  is  a  plot  of  U  for  three  different  field 
strengths.  From  this  plot,  it  is  found  that  U  is  clamped 
by  the  field  of  the  standing  wave,  and  that  as  the  amplitude 
of  the  wave  increases,  XI  is  clamped  down  to  zero.  It  can 
also  be  seen  that  at  the  nodes  of  the  standing  wave,  \xi\ 
relaxes  to  a  larger  value.  Upon  closer  observation,  the 
tail  immediately  following  the  standing  wave  node  is  found 
to  be  larger  than  the  one  preceding  it.  Figures  (2)  and 
(3)  show  that  the  radiation  force  is  negative  preceding  each 
node,  and  positive  immediately  following  them. 


5  6  7  8  9 

time  (*  I0'^5cc) 


Figure  5.  Bloch  variable  I*,  and  force  (dynes)  for 

IA  around  a  node  in  the  standing  wave. 
(Shaded  area  is  the  difference  between  the 
areas  under  the  positive  and  negative 
force  curves) . 


time  (x  IO'S5ec.) 


Figure  6.  Bloch  variable  ti  and  force  (dynes)  for 

Ho3  5 A  around  a  node  in  the  standing  wave. 
(Shaded  area  is  the  difference  between  the 
areas  under  the  positive  and  negative 
force  curves) . 


The  product  of  *U  and  s»n  kvt  is  on  the  order  of  unity/ 
so  the  magnitude  of  the  force  is  approximately 


F|  ~  i-Kkn.. 


The  direction  of  this  force  depends  on  the  portion  of 
the  standing  wave  that  the  atom  is  in.  This  force  acts  for 
the  distance  of  approximately  . 

If  the  atom  is  moviny  with  a  constant  velocity  V  , 
then  the  length  of  time  that  the  force  will  be  acting  on  the 
atom  is  r  =  . 

The  change  in  momentum  of  the  atom  will  then  be 

kp  «  mSv  =  Ft  (49) 

or 


(50) 


The  percentage  change  in  velocity  is  thus  given  by 

w  _  ir-tSilo 

or,  by  using  approximate  values, 

fev  _  OOPo*)  _  J_ 

V  (10  10* 


(51) 


(52) 


For  the  velocities  considered  in  this  study 


IO*  to 


|0  tm/sec  ),  the  change  in  velocity  over  1/4  period  of  the 
standing  wave  is  on  the  order  of  .01%.  Thus,  the  assumption 
of  a  constant  velocity  is  a  fairly  good  approximation.  If, 


however,  the  atomic  velocity  is  of  the  order  of  10  or 
less,  or  the  field  strength  is  of  the  order  of  £1,-  10* 
to  |0°  ,  then  the  constant  velocity  assumption  is  of  a 

questionable  value. 


Variable  Velocity  Program 

To  get  an  accurate  picture  of  the  atomic  motion  in  a 
standing  wave,  the  velocity  of  the  atom  cannot  be  assumed  to 
be  constant.  It  will  vary  according  to  the  strength  of  the 
light  pressure  the  atom  encounters  and  the  length  of  time 
the  force  is  applied. 

The  computer  program  differs  from  the  constant  velocity 
Ehren test- Bloch  program  in  that  the  complete  Ehren f est-Bloch 
equation  set  is  used  to  calculate  the  motion  of  the  atom. 

The  instantaneous  force  is  calculated  at  each  time  step, 
then  used  to  calculate  the  change  in  atomic  velocity.  The 

.  -23 

atomic  mass  was  taken  to  be  rn  —  *  10  g  which  represents 

an  atom  about  the  weight  of  a  sodium  atom.  Again,  the 
A-coef f ic ient  and  the  resonant  transition  frequency  were 
chosen  as  convenient  values  only,  and  do  not  represent 
values  for  real  atoms. 

The  average  light  pressure  force  on  the  atom  can  be 
calculated  by  the  average  change  in  the  velocity  over  one 
period  of  the  standing  wave. 

—  kv 

F  =  m  (53) 


Several  points  of  the  average  force  were  calculated  this  way 
and  are  plotted  as  triangles  against  the  calculated  curve  in 


Fast  Fourier  Transform 


A  256-point  Fast  Fourier  Transform  (FFT)  program  was 
written  in  Fortran  77  based  on  a  program  for  a  microcomputer 
[12],  This  program  was  used  to  decompose  the  functions  for 
the  Bloch  variable  'Ll  and  the  radiation  force  from  the 
constant  velocity  program,  into  their  constituent  harmonic 
frequencies.  The  time  increment  was  adjusted  to  allow  256 
time  steps  per  period  of  the  standing  wave,  and  the  output 
of  the  FFT  program  was  a  descrete  harmonic  frequency 
spectrum  of  the  variable  being  studied.  This  program  is 
listed  in  Appendix  B. 

Figure  (10)  shows  the  frequency  spectrum  for  the  Bloch 
variable  U.  for  several  field  strengths.  As  expected  from 
studying  the  weak  field  case,  the  contribution  by  the  higher 
harmonic  terms  is  dependent  on  the  intensity  of  the  field, 
with  the  higher  harmonic  terms  becoming  negligible  for 
relatively  weak  fields. 

The  radiation  force  frequency  spectrum  plotted  in 
Figure  (11)  also  shows  the  dependence  of  the  higher  harmonic 
terms  on  the  field  intensity.  The  force  is  seen  to  contain 
a  set  of  even  harmonic  terms  in  accordance  with  the  results 
of  Minogin  and  Serimaa  in  equations  (2).  The  use  of  this 
FFT  program  is  limited,  however,  because  of  the  inability  to 
match  exactly  256  steps  in  one  period.  This  caused  the 
harmonic  terms  to  bleed  over  onto  the  adjacent  harmonics. 

To  correct  for  this,  it  would  be  necessary  to  run  a  512 
point  or  a  1024  point  FFT  program  on  the  same  interval. 


Figure  11. 


Descrete  harmonic  frequency  spectrum 
for  the  radiation  force  calculated 
at:  (a)  _Q.„=|A  ,  (b)il„=5A  , 

and  (c)  SI.  =  10  A. 


Conclusions  and  Recommendations 


This  study  of  atomic  motion  in  a  standing  wave  was 
conducted  in  two  parts:  The  first,  an  analytical  approach 
to  solving  the  optical  Bloch  equations,  and  the  second,  a 
numerical  approach  to  solve  the  complete  Ehrenfest-Bloch 
equations-  The  analytical  approach  was  unsuccessful  in 
obtaining  a  complete  solution  to  the  Bloch  equations,  but 
several  recursion  formulas  were  obtained  that  can  be 
approximated  and  solved  for  specific  cases  such  as  the  weak 
field  case.  The  numerical  approach  was  successful  in 
obtaining  instantaneous  values  for  the  Bloch  variables  and 
the  light  pressure  force.  Several  examples  of  calculations 
done  numerically  were  shown  to  illustrate  the  use  of  the 
computer  programs. 

From  the  numerical  calculations,  it  was  shown  that  the 
Ehrenfest-ULoch  equations  indeed  represent  the  same  elements 
of  resonance  radiation  pressure  theory  as  the  theories 
developed  by  Minogin,  Letokhov,  and  Serimaa.  This  is 
because  both  theories  are  developed  along  the  same  lines 
from  the  dipole  representation  of  the  interaction  between 
the  atom  and  the  electromagnetic  wave.  The  difference  is 
that  Cook's  theory  goes  a  step  farther  to  use  the  Bloch 
vector  representation  of  the  atomic  density  matrix  to  make 
the  theory  somewhat  simpler  and  easier  to  visualize.  The 
resulting  Ehrenfest-Bloch  theorem  allows  for  a  quick 
interpretation  of  the  internal  3tate  of  the  atom  during  the 
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interaction  with  the  light  wave. 

The  Fourier  series  solution  obtained  by  Minogin  and 
Serimaa  [11]  can  probably  be  obtained  from  the  optical  Bloch 
equations,  but  because  of  the  complexity  of  the  convergent 
continued  fraction  coefficients,  its  usefulness  beyond  some 
reasonable  approximation  is  questionable.  With  an  extremely 
simple  computer  program,  the  Ehrenfest-bloch  equations  can 
be  solved  numerically  to  get  the  same  results  as  shown  by 
comparing  Figures  (1)  and  (7). 

Another  difficulty  with  the  solution  developed  by 
Minogin  and  Serimaa  lies  in  the  assumption  of  a  constant 
atomic  velocity.  It  was  shown  that  for  low  to  moderate 
field  strengths,  this  assumption  was  valid,  but  as  the  field 
strength  increased,  the  assumption  of  constant  velocity 
became  a  questionable  approximation.  In  contrast,  the 
Ehrenf est-Bloch  equations  can  be  programmed  to  allow  for  a 
variable  atomic  velocity,  a  feature  that  would  be  necessary 
for  very  slow  atoms  and  very  intense  fields. 

One  interesting  feature  which  was  found  by  Minogin  and 
Serimaa  and  explained  by  Kyrola  and  Stenholm  was  the 
existence  of  force  variations  including  a  force  reversal  for 
relatively  slow  atomic  velocities.  The  explaination  of  an 
multi-photon  or  doppleron  resonance  with  the  two  counter- 
propagating  light  waves  in  the  standing  wave  makes  good 
physical  sense.  But  by  investigating  the  phenomenon  using 
the  Bloch  vector  representation,  it  may  be  possible  to 
develop  a  theory  of  resonance  involving  the  classical  dipole 


moment  of  the  atom  and  the  classical  electromagnetic  wave. 
Such  a  resonance  could  very  likely  include  some  form  of  the 
Rabi  flopping  frequency  which  describes  the  absorption  and 
subsequent  emission  of  energy  from  the  wave. 

Thus,  this  thesis  leads  to  several  recommendations  to 
further  the  study  of  resonance  radiation.  First,  the  Bloch 
representation  of  the  atom  can  be  studied  in  the  standing 
wave  interaction  to  try  to  understand  and  explain  the 
resonance  that  causes  reversals  in  the  average  light 
pressure  force.  Second,  the  Ehrenf es t-Bloch  equations  can 
be  programmed  to  represent  an  actual  experimental  setup  and 
the  determine  the  results  that  are  expected.  Third,  the 
analytical  solution  can  be  pursued  to  obtain  a  solution 
similar  to  equations  (2),  but,  as  stated  before,  the  use  of 
such  a  solution  may  have  limited  use  except  to  verify  the 
theoretical  development. 

An  area  that  deserves  more  investigation  is  when  the 
time-averaged  light  pressure  force  changes  sign  for 
relatively  small  atomic  velocities.  Because  the  force  goes 
to  zero  at  three  atomic  velocities,  v=  O  and  va  *  V,  , 
the  effect  on  a  distribution  of  atomic  velocities  would  be 
to  separate  the  atoms  into  three  groups,  one  that  has  zero 
velocity  and  two  that  are  moving  in  opposite  directions  with 
the  same  speed.  Such  an  experiment  may  be  accomplished 
using  a  slow  atomic  beam  passing  perpendicularly  through  a 
large  region  with  a  standing  laser  wave.  The  transverse 
thermal  velocities  will  be  damped  out  and  the  resulting  beam 


profiie  whould  show  the  three  transverse  velocity  groups. 
Another  experiment  could  involve  the  use  of  a  cavity 
resonator  in  which  a  distribution  of  atoms  would  be 
irradiated,  then  the  spectrum  immediately  observed  with  a 
probe  laser  to  determine  the  velocity  distribution  of  the 
atoms. 

Through  a  judicious  use  of  resonant  radiation  pressure 
it  is  possible  to  manipulate  matter  on  the  atomic  level. 

Such  experiments  can  be  applied  to  test  theoretical 
calculations  of  light  pressure  and  may  eventually  lead  to 
the  successful  cooling  and  trapping  of  neutral  atoms. 


APPENDIX  A:  The  Optical  Bloch  Equations 

The  optical  Bloch  equations  (A-l)  form  a  system  of 
three  equations  which  completely  describe  the  internal 
motion  of  a  two-level  atom  interacting  with  light  in  the 
Rotating  Wave  Approxima tion . 

u  -  -lA-u  +  ^'ir 
ir  =  -  ^Aar  •+•  5L-ur 


'W  -  -SLnr  -  A(w  +  \) 


( A-la) 
( A-lb) 
( A-lc) 


Y  =  (ti.ry)  is  the  Bloch  vector  of  the  atom,  A  is 
the  Einstein  A-coef f icient  for  spontaneous  emission,  SL  is 
the  Rabi  flopping  frequency,  and  &  is  the  detuning  of 
the  frequency  of  the  incident  light  from  the  resonant 
transition  frequency  of  the  atom. 

The  Bloch  equations  are  derived  from  a  representation 
of  the  equation  of  motion  of  the  density  matrix,  ,  of 

the  atom.  The  density  matrix  is  defined  by 


(3  =  z  p„  now 


(A-2) 


where  X  is  the  wavefunction  of  the  state  fX  ,  and 
is  the  probability  of  the  atom  being  in  that  state.  The 
equation  of  motion  is  given  by  the  Heisenberg  equation: 


If-  as [M] 


at 


( A-3a ) 
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or  in  matrix  form, 


f„w=  k  ~  P„kHKm)  (ft-: 

Each  state  has  its  corresponding  eigenvalue  such  that 

H  —  EK£-m  (A- 

1  nm  h  nr*\  ' 

Equation  (A-3b)  thus  becomes 

Km  “  ^  “  Pnk^k^km)  (A 

Km  ~  7* 

or 

f_  = 

Now,  the  Bohr  transition  frequency  between  two  states  is 


thus  equation  ( A-7 )  becomes 


C  —  -fco,  p  (A-9) 

\  r\m  ft"'\  nm 

The  solution  of  this  equation  is 

*Jt)  *  W  “-10> 

In  the  case  of  n  =  m  ,  then  <Jnyn,=  0  ,  and  the  probability 
of  being  in  state  n  would  be  time  independent.  This  is 
corrected  by  adding  the  effects  of  spontaneous  emission. 

For  diagonal  matrix  elements  (  n  ■=.  m  )  ,  then 


equation  ( A-9 )  becomes 


whe  re 


]C  ^na^nr\  (A 

p  >n  '  ^<11  » 

e  A„m-  ^nw./^'Kc3  .  The  first  term  in 


(  A— 1 1  ) 


equation  (A-ll)  represents  the  increase  in  probability  of 
being  in  state  n  due  to  the  decay  from  states  p  above 
state  n  ,  and  the  second  term  represents  the  decrease  due 
to  losses  to  lower  states 

For  off  diagonal  matrix  elements  (n  ^  rn )  ,  equation 
(A-9)  becomes 

A„,  +  Z  (A*12) 

'  »  ' 

Now,  it  the  atom  interacts  with  an  external  electric 


field,  the  Hamiltonian  can  be  represented  by 

A  A  A  , 

H  =  H  +  H 


(A-13) 


A  I  A 


where  H  -  —  yd  •  E-C't') 


.  If 


£(t)  =■  £  E  cosot  i 


is  the 


electric  field  vector  for  a  light  wave  of  frequency  Co  , 
then 


~  '  A _ 

H  -  toi  ot 


(A-14  ) 


A  -  ' 

where  pt-  fe 

is  the  dipole  moment  operator  of  the  atom 
Equation  (A-9)  becomes 


(A-15) 


With  the  spontaneous  emission  terms,  this  becomes 

l 


^  nr*  ~  ~  ^^V') 

^  ^pnPpP  H  ^n^Cni 

p>n  1  1 


for  r\  —  m  ;  and 

P  =  "  +  ^  ^  ^  ?‘lm 

~"  11  (  £«  ^  n  *  £L^"0?r 


,  r\r»\ 


■  ||«K  V  <|4|n 

for  r\  rr\ 

For  a  two-level  atom,  the  density  matrix  is  a 
matrix 

f„  P, 

?*  e. 


f. 


I* 


on* 


ix 


and  the  dipole  moment  operator  matrix  is 


<r>  - 


o  f 


r- 


o 


where  J-*-  is  a  real  number.  Equations  (A-16a)  and 
become 

P„=  rir("rEco:>^t  -*■  rEo,stot 

4  j^Ecoswtp^  -  A„pw 

?,*  -  -  “*>.« + ih (')*E wt pi, +  rE °* ^ p4  -  i  ^zv p.o. 

*  # 

e  =  p 

'  Z»  '  I 


0  11 
IX 


l A-16a ) 

(A-16b) 

hermitian 

( A— 17  ) 

(A-18) 

( A-16b) 

( A-19a) 

( A-19b) 

( A-19  c ) 

( A-19d) 
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p„=  +  Ae„ 

p„=  ex)  -  Ae„ 

Prt  =  io„  4  iSlcoSW*  fe.-p„) 

f*.  =  e,; 


( A-20a ) 
( A-20b ) 
( A-20c ) 
( A-20d) 


where  Si-  ^e/-^  is  the  Rabi  flopping  frequency  and  A®  Ax»  . 
In  these  equations,  —  tOt,  "=■  —  iOilm  is  the  resonant 

transition  frequency  of  the  two-level  atom. 

Now,  the  density  matrix  of  the  atom  can  be  rewritten  in 
terms  of  a  new  matrix  0"nm  that  is  a  slowly  varying  matrix 
in  the  near-resonance  case  in  which  the  dominant  term  (  co 
near  C00  )  has  been  removed. 


-  <r.  ^ 

Equations  (A-20)  then  become 


(A-21a) 
( A-21b ) 
( A-21c) 


%  = 

-iSlam#  ( trie“A,t-  <j^  e~l“>  ) 

*-  A  cr^ 

( A-22a ) 

<Slx~ 

iSl  c«^(o;xeiwt'a-2<e^) 

-  A  <rtl. 

( A-22b) 

=  iSlco»u>t(<r^-<rll)t.' ,ut 

—  C  (  u3  - 

( A-22c) 

o1  * 

i-r 

( A-22d) 

4’ 

s  i 


co%  ot  e  -  2.  \  +  e 


x  e 


IJ  ( A— 23 ) 


In  the  Rotating  Wave  Approximation  (RWA),  the  terms  which 


vary  at 

twi 

Lee  the 

dom 

inant  f 

vary  ra 

pidly  and  . 

average  to 

t  co+" 

1 

Cos 

cot 

e 

-JV 

~L 

Therefore , 

equa  t  i 

ons 

( A-2  2  ) 

% 

=  - 

±SL  c 

(<ni 

-«!.) 

— 

-iSL  i 

1 

zs 

ksi  i( 

°zv 

02, 

— 

Now,  the  Bloch  vector  is  defined  by 

'U  =  -t-  -  2  (Re  [o'.t] 

'V  =  -  t  (cr;v  t  (fj  ~  2ci3«i  [ 

'ViT  =.  a;.-  <ru  ^  P,  -  P. 


( A— 24  ) 


( A-25a) 
( A— 25b ) 
( A-25c ) 
( A-25d ) 


( A-2t>a ) 
( A-26b) 
( A-26c ) 


where  the  variable  'W  is  the  inversion  between  the  two 
states.  In  the  case  of  the  two- level  atom,  P  •*-  —  | 


k 


Subtracting  (A-25a)  from  (A-25b)  provides 


i 

*],)  =  &  ‘  -  A  (o"xx  + 

l  A-27  ) 

(<f„-  iTu)  =  -Sl[-i(<rtt-<rJ]  -A[(<V'i;.)*-(<rll,0] 

( A-28 ) 

•>  . 

1” 

i 

L- 

Adding  (a-25c)  to  (A-25d)  provides 

,  4 

(^iz ~~  ~  ‘xA<^zl  ~  i  -V  c  A  <3"^, 

( A— 29 ) 

■v:1 

P 

r/ 

>. 

+A[-i(<r  -O] 

( A-30  ) 

,  % 

\ 

K 

r 

And  subtracting  (A-25d)  from  (A-25c)  provides 

p 

(oU~  <*V)  -  ~  -  4:Ac(71  +  iAcrxl 

*„ 

{ 

+  ^.i) 

( A-31 ) 

| 

s 

-USx-  «v)  =  -  -sA[-i(<j'.x-«v.')] 

r- 

7~i 

■ 

+  S}.(a-m.-<rll') 

l A-32) 

V 

B 

-J 

And  finally,  with  the  appropriate  substitutions,  the 

familiar  form  of  the  optical  Bloch  equations  is  derived 

y/y.\ 

V  1 

• 

*  ' 

=  -  k  h'u  +-  Mr 

( A-33a) 

- 

/« 

» 

"  ' 

or  —  ~  A.'U  -  i  A  v  +  51 'ur- 

( A-33b) 

!•;-  .• 

• 

ad"  —  —  Sl'lr  -  A (ov~  +-  l') 

( A-33c) 

r 

1 

The  Bloch  vector  provides  a  useful  description 

of  the 

-> 

• 

rr 

two-level  atom  in  the  RWA.  IX  and  'V~  are  components 

r- 

•*-V. 

. *  ' 

• 

^  Sj 

- 

4 

* 

APPENDIX  B:  Computer  Source  Code 


C  BLOCH  SOLVER 

C 

C  THIS  PROGRAM  CALCULATES  THE  BLOCH  VECTOR  COMPONENTS  AND 
C  THE  INSTANTANEOUS  LIGHT  PRESSURE  FORCE  ON  AN  ATOM  USING 
C  THE  EHKENFEST— BLOCH  EQUATIONS  WITH  A  CONSTANT  VELOCITY 
C  ASSUMPTION. 

C 

CHARACTER *8  FILEl 
REAL  A/C/H/K/T/ DEL/FO/Fl/ OM 
REAL  UO / VO ,  WO , FRC , V 
COMMON  K  ,  V  , DEL /OM 
C 

C=2.99792458E10 
F0=1 . 0E14 
DEL=— 1.0E8 
Fl=FO+DEL 
K=F1/C 
C 

C  THE  ATOM  STARTS  IN  THE  GROUND  STATE  (0,0,-l) 

C 

l’=0 . 0 
U0=0.0 
V0=0.0 
W0=-1.0 
C 

WRITE (6, 5) 

5  FORMAT ( / '  ENTER  OUTPUT  DATAFILE  NAME’) 

READ ( 6 / 10 )  FILEl 
10  FORMAT (A6) 

C 

OPEN ( UN IT =10 /FILE=FILEl ,STATUS= ' NEW' , FORM= ' FORMATTED' ) 
C 

WRITE (6*15) 

15  FORMAT ( /// '  INPUT  OMEGA') 

READ (6, 20)  OM 
20  FORMAT (E6. 4) 

WRITE (6/25 ) 

25  FORMAT ('  INPUT  ATOMIC  VELOCITY’) 

READ (6, 20)  V 
C 


C  STORE  THE  RUN  INFORMATION  IN  FILE  10. 

C 

40  WRITE (10 ,50)  OM, V , FO, Fl , DEL 

50  FORMAT (////'  OMEGA  =  ',E10.5/ 

+  '  ATOMIC  VELOCITY  =  • ,E10.5/'  TRANSITION  ', 

+  'FREQUENCY  =  ',E10.5/'  LASER  FREQUENCY  =  ', 

+  E10.5/’  DELTA  =  ' , Ell . 5/ / ) 


WRITE (10 ,60  ) 

FORMAT (6X, 1  TIME'  ,9X, ' U( t) '  ,  10X, ' V( t) '  ,10X, ' W( t) '  ,10X, 

+  '  FORCE  '  /  ' - '  , 

+  < - •  ) 

WRITE (10, 110)  T , UO , VO , WO , FRC 

CALCULATE  THE  BLOCH  VECTOR  UNTIL  THE  STOP  TIME  AND  STORE 
THE  VALUES  IN  FILE  10. 

CALL  BLOCH (T , UO , VO , WO , FRC ) 

WRITE ( 10 , 110  )  T,UO,VO,WO,FRC 
i  FORMATC  '  ,E10.3,4(2X,E12.5)) 

IF(T.LE.5.0E— 6)  GOTO  100 

STOP 

END 


SUBROUTINE  BLOCH ( T , UO , VO , WO , FRC ) 
REAL  T , UO , VO , WO , FRC 
REAL  K , V , DEL , OM 
COMMON  K , V , DEL ,0M 

A=1 . 0E8 

H=1 .OE-11 

HBAR=1 .054592E-27 

DO  200  1=1,500 

OMEGA=OM*COS ( K*V*T ) 

OPTICAL  BLOCH  EQUATIONS 

0U=— . 5*A*U0+DEL*V0 

DV=— DEL*UO— . 5*A*VO+OMEGA* WO 

DW=-OMEGA* VO— A* WO— A 


CALCULATE  THE  BLOCH  VECTOR  AFTER  THE  TIME  STEP  H. 


T=T+H 

UO=DU*H+UO 
VO=DV*H+VO 
W0=DW*  H+WO 

)  CONTINUE 

EQUATION  OF  MOTION  FROM  EHRENFEST-BLOCH  EQUATIONS 
FRC=- . 5*K*HBAR*U0*OM*SIN ( K*V*T ) 


RETURN 


v.v.v.’y.vj  w:w>,9  \w  j »» *  gm  *  ■  n 


Variable  Velocity  Program 


C  EHRENFEST-BLOCH  SOLVER 
C 

C  THIS  PROGRAM  CALCULATES  THE  BLOCH  VECTOR  COMPONENTS  AND 
C  THE  INSTANTANEOUS  LIGHT  PRESSURE  FORCE  ON  AN  ATOM  USING 
C  THE  EHRENFEST-BLOCH  EQUATIONS  WITH  THE  VELOCITY  ALLOWED 
C  TO  VARY. 

C 

CHAKACTER*8  FILEl 

REAL  A,C,H,K,T, DEL, FO,Fl,OM, MASS 
REAL  UO, VO,WO,FRC, X, V 
COMMON  K, DEL, OM, MASS 
C 

MASS =4 . OE-23 
C=2 . 99792458E10 
F0=1 .0E14 
DEL=-1 . 0E8 
F1=F0+DEL 
K=F 1/ C 
C 

C  THE  ATOM  STARTS  IN  THE  GROUND  STATE  (0,0,-l) 

C 

T=0 . 0 
X=0 . 0 
U0=0.0 
V0=0.0 
W0=~1 .0 

c 

WRITE (6, 5) 

5  FORMAT ( / 1  ENTER  OUTPUT  DATAFILE  NAME') 

READ (6,10)  FILEl 
10  FORMAT (A6) 

C 

OPEN ( UN IT =10 ,FILE=FILEl ,STATUS= ' NEW' ,FORM=' FORMATTED' ) 
C 

WRITE (6,15) 

15  FORMAT (/’  INPUT  OMEGA') 

READ (6, 20)  OM 
20  FORMAT (E6. 4) 

WRITE (6, 25 ) 

25  FORMAT (’  INPUT  INITIAL  ATOMIC  VELOCITY') 

READ (6, 20)  V 
C 

C  STORE  THE  RUN  INFORMATION  IN  FILE  10. 

C 

40  WRITE (10, 50)  OM , V , FO , Fl , DEL 

50  FORMAT ( / // / '  OMEGA  =  ' ,El0.5/ 

+  '  INITIAL  ATOMIC  VELOCITY  =  '  ,El0.5/’  TRANSITION  '  , 

+  'FREQUENCY  =  *,E10.5/'  LASER  FREQUENCY  =  ', 

+  E10.5/'  DELTA  =  ' ,E11 . 5//) 

C 
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WRITE (10, 60) 

FORMAT (6X, 'TIME'  ,9X, ' U( t) '  , 10X, ' V( t) '  , 10X, ' W( t) '  ,9X 
' FORCE*  , 9X , ' X( t) '  ,9X,  ' VEL(  t) ' / * - 


WRITE (10,110)  T , UO , VO , WO , FRC , X , V 

CALCULATE  THE  BLOCH  VECTOR  UNTIL  THE  STOP  TIME  AND  STORE 
THE  VALUES  IN  FILE  10. 

CALL  EH REN (T,UO,VO,WO,X,V,FRC) 

WRITE (10, 110)  T,UO, VO, WO, FRC, X, V 
FORMAT ('  ' ,E10.3, 5(2X,E12.5)  ,2X,E15.8) 

IF(T.LE. 5.0E-6 )  GOTO  100 

STOP 

END 


SUBROUTINE  EHUEN ( T , UO , VO , WO , X , V , FRC  ) 

REAL  T, UO, VO, WO, FRC, X, MASS 
REAL  K , V , DEL ,OM 
COMMON  K, DEL, OM, MASS 

A=1 . 0E8 

H=1 .OE-11 

HBAR=1 . 05459 2 E— 27 

DO  200  1=1,500 

OMEGA=OM*COS(K*X) 

OPTICAL  BLOCH  EQUATIONS 

DU=— . 5*A*U0+DEL*V0 
DV=-DEL*UO— . 5*A* VO+OMEGA* WO 
DW=— OMEGA* VO— A* WO— A 

CALCULATE  THE  BLOCH  VECTOR  AFTER  THE  TIME  STEP  H. 
T=T+H 

UO=DU*H+UO 

VO=DV*H+VO 

WO=DW*H+WO 

X=X+V*H 

EQUATION  OF  MOTION  FROM  EHRENFEST-BLOCH  EQUATIONS 

FRC=-. 5*HBAR*K*UO*OM*SIN (K*X ) 

V=V+FRC*H/ MASS 

)  CONTINUE 

RETURN 

END 
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PROGRAM  FASTRAN 


A  256-point  FAST  FOURIER  TRANSFORM  program  adapted  from 
a  program  which  appeared  in  BYTE  Magazine  in  December 
1978. 


CHARACTER*8  FILE1 , FILE2 
REAL  XI (256 ) ,X2 (256  ) 

REAL  UO ( 256 ) , VO (256 ) ,W0(256 ) ,FRC(256  ) 

INTEGER  X , Y / U 
COMMON  L / N 
N  =  25b 
l=b 

PI=3. 141592654 

INPUT  TIME-DEPENDENT  FUNCTION 
WRITE (6, 5) 

FORMAT (/'  ENTER  INPUT  DATAFILE  NAME’) 

READ (6,6)  FILEl 
FORMAT ( A6  ) 

WRITE (6, 7) 

FORMAT ( / '  ENTER  OUTPUT  DATAFILE  NAME*) 

READ (6, 6)  FILE2 

WRITE (6, 10  ) 

FORMAT ( / '  ENTER  STARTING  LINE  NUMBER')  * 

READ (6,15  )  INDO 
FORMAT ( 13  ) 

OPEN(UNIT=10,FILE=FILE1 , STATUS= ' OLD ' , FORM= ’ FORMATTED ' ) 
OPEN ( UNIT =11 ,FILE=FILE2 ,STATUS=' NEW' ,FORM=’ FORMATTED' ) 

DO  30  1=1 , INDO— 1 
READ (10, 20)  XXX 
FORMAT (Al) 

CONTINUE 

DO  50  1=0,256 

READ (10,40)  UO ( I ) , VO ( I ) , WO ( I ) , FRC ( I ) 

FORMAT (T14 , E 12 . 5 , T28 , S 12 . 5 ,T42 , E12 . 5 ,T56 , El2 . 5  ) 

CONTINUE 

WRITE (6, 60  ) 

FORMAT ( / '  WHICH  VARIABLE  DO  YOU  WANT?  (1  -  4)') 

READ (6, 70)  IV  AR 
FORMAT (II ) 

GOTO (75,80,85,90)  IVAR 
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DO  76  1=0, N 
X1(I)=U0(I) 

CONTINUE 
GOTO  99 

DO  81  I =0  ,  N 
XI  (I)=V0  (I) 

CONTINUE 
GOTO  99 

DO  86  I  =0  ,  N 
XI (I)=W0  (I) 

CONTINUE 
GOTO  99 

DO  91  1=0, N 

Xl(I)=FRC(I)*.lEl7 

CONTINUE 

CONTINUE 

SCALE  INPUT  TIME  FUNCTION 

DO  100  1=0, N-l 
XI  (I)=X1 (I)/N 
)  CONTINUE 

FFT  IN-PLACE  ALGORITHM 

I l=N/2 
12=1 

V=2*PI/N 

DO  400  1=1, L 

13  =  0 

14=11 

DO  300  J=1 , 12 
X= I  NT ( I  3/11 ) 

CALL  SCRAM (X,Y) 

I  5=  Y 

Zl=COS (V* I  5 ) 

Z2=-SIN (V* 1 5 ) 


DO  200  K=I 3 , 14-1 
Al=Xl (K) 

A2=X2(K) 

Bl=Zl*Xl (K+I1)-Z2*X2 (K+Il) 
B2=Z2*Xl (K+Il ) +Z1*X2 (K+Il) 
XI  (K ) =A 1+Bl 
X2  (K) =A2+B2 
XI  (K  +  Il) =A1-B1 
X2  ( K+ 1 1 ) =A2-B2 


13=13+2*11 
14=14+2*11 
300  CONTINUE 
C 

11=11/2 

12=2*12 

400  CONTINUE 
C 

C  OUTPUT  TABLE  OF  VALUES 
C 

U=0 

C 

IF( IVAR.EQ. 1  )  WRITE (11/410) 

IP( IVAR.EQ.2)  WRITE(11 ,420 ) 

IF( IVAR.EQ. 3)  WRITE (11/ 430 ) 

IF( IVAR.EQ. 4 )  WRITE (11 ,440 ) 

410  FORMAT (/// '  FFT  for  U(t)’) 

420  FORMAT (///'  FFT  for  V(t)') 

430  FORMAT (///'  FFT  for  W( t) ' ) 

440  FORMAT (///'  FFT  for  FRC(t)') 

C 

WRITE (11 ,450 ) 

450  FORMAT ( / '  HARMONIC  REAL  IMAGINARY  MAGNITUDE') 

WRITE (11, 460) 

460  FORMAT  ('  - ’) 

C 

500  X=U 

CALL  SCRAM (X,Y) 

X3=SQRT(X1 (Y)*Xl(Y)+X2(Y)*X2(Y) ) 

WRITE (11, 550 )  U,X1(Y) ,X2(Y) ,X3 
550  FORMAT ('  '  , 14 , 3 (IX, F10 .6  ) ) 

U=U+1 

IF ( U . LE . N/2 )  GOTO  500 
C 

WRITE (6, 600 ) 

600  FORMAT (/' CALCULATE  FFT  FOR  ANOTHER  VARIABLE?') 

READ ( 6 , 650 )  YES 
650  FORMAT (Al) 

IF(YES.EQ. ' Y' )  GOTO  55 
C 

CLOSE (10) 

CLOSE ( 11  ) 

C 

STO  P 
END 
C 
C 

C  SUBROUTINE  SCRAMBLER 
C 

SUBROUTINE  SCRAM(X,Y) 

INTEGER  X , Y 
COMMON  L , N 
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